Вестник ДГТУ. 2011. Т. 11, № 6(57) 








УДК 539.3 


ТРЕХМЕРНАЯ КОНТАКТНАЯ ЗАДАЧА О ВЗАИМОДЕЙСТВИИ УПРУГОГО СЛОЯ 
С ДВУМЯ ШТАМПАМИ ПРИ УЧЕТЕ ТРЕНИЯ” 


Д.А. ПОЖАРСКИЙ, А.А. МОЛЧАНОВ 
(Донской государственный технический университет) 


Изучены пространственные контактные задачи для упругого слоя конечной толщины, в одну грань которого 
симметрично вдавливаются два одинаковых жестких эллиптических штампа с учетом трения при разных ти- 
пах граничных условий на друг ОЙ Г рани. Задачи сведены к интеп ральным уравнениям относительно контакт- 
ного давления, которые решены методом Галанова. 
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Введение. Исследованы трехмерные контактные задачи теории упругости при учете трения о 
взаимодействии слоя с двумя симметричными штампами, расположенными на одной его грани. 
Другая грань слоя находится в условиях жесткой или скользящей заделки. Штампы имеют форму 
эллиптических параболоидов, начинают удаляться друг от друга или сближаться. Области кон- 
такта неизвестны. Ранее аналогичные задачи с трением рассматривались для случая одного 
штампа на слое [1 - 4] и на полосе [5]. В работе [6] рассматривается пространственная задача о 
контакте с упругим слоем системы двух симметричных эллиптических штампов с плоской подош- 
вой (асимптотический метод решения). 


Постановка задачи. В декартовых координатах рассмотрим слой 1. у е (—, о), Е [0, #]} тол- 


щиной А, нижняя грань которого 2 = 0 находится в жесткой или скользящей заделке (задачи А и Б 
соответственно). Упругий материал слоя имеет коэффициент Пуассона у и модуль сдвига С. Верх- 
няя грань слоя 2 = А взаимодействует с двумя одинаковыми штампами. Формы основания штампов 
имеют вид эллиптических параболоидов и описываются функциями 
9+(жу) = (же /(281)+У 28»), В > К. (1) 

Между поверхностью слоя и штампов действуют силы кулоновского трения с коэффициен- 
том трения и. Штампы начинают достаточно медленно двигаться вдоль оси х так, что задачи 
симметричны относительно оси у. Силы трения направлены против движения. При в > 0 штампы 
удаляются друг от друга, а при в < 0 — начинают сближаться. К штампам приложены симметрич- 
ные касательные силы Т, нормальные силы Р. Пусть осадка штампов равна 5, а перекос отсутст- 
вует. Симметричные по у области контакта ©. неизвестны (область ©. при х < 0). 

При известных величинах С, у, ц, А, с, В, Е и би заданной функции Кг,2) требуется оп- 
ределить контактные давления с.(х, у, А) = —-В(х,у), (х,У)ЕО = ОО, а также сами области контак- 
та О.. Затем можно найти, например, величину Риз условия равновесия штампа 


РЕ Гас, уахау . (2) 
о, 


Из аналогичных интегральных условий можно найти величину Т, а также плечи приложе- 
ния сил Ри Т. 
Решение задачи. Предположим, что области контакта априори содержатся в прямоугольниках 
5 = {+ <а, У<Ь}, Б>а, с>а. 

Для вывода интегральных уравнений (ИУ) контактных задач А, Б используется интеграль- 
ное преобразование Фурье и закон Кулона. В результате придем к ИУ на двух участках контакта, 
которое после введения безразмерных обозначений 


р с а с 6 Ь Ь 
- ПЬ в ‚б= ‚ А= ‚ В 2 (3) 
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* Работа поддержана грантом РФФИ 09-01-00004. 
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4(х, у) т 30. 
(ху) = ‚ Р'= ‚ 0= ‚ О'<>0, О. <>О0,, 5, <>5, 4 
об ВЫ ЕТ. (9) 
можно записать в виде (штрихи далее опускаем): 
1 х-ё 1 х-& ут м 
Г Г т ТЕ ЕТ ват = Л(х,у), (х,у)ЕО, (5) 
ЛО, у) =6- А(х+с)* - Ву’, (х,У)ЕС,, (6) 
К, =[(< +8)" +()-п)?]°, а=на-2у)/(2-2У), (7) 


Т.о = [о (и) - Пи? +1? )аи+ ть ‚(и)- ПЛ (и + т? )аи. (8) 
Здесь (и) — И Бесселя. В ядре ИУ (4) выделена главная часть. 
Для задачи А (жесткая заделка) 
2к $12и-4и 
Г = ‚ КЕ=З-—4у, 9 
"() 2к си2и +4и? +1+к? 9 


2к сВ2и- 4(1- 2%) и? -2к 











Г. = 10 
2 (и) 2к с62и+4и? +1+к? (о 
а для задачи Б (скользящая заделка) 
сви —1 $62и - 21-2) 'и 
(и) =———, [,(и)= . 11 
‚() $62и+2и 2 (и) $62и + 2и с. 


Введенный в (3) безразмерный параметр ^, характеризует относительную толщину упруго- 
го слоя с учетом симметрии задач а(-х.у) = а(<ху). Тогда уравнение (5) сводится к ИУ на одном 
участке контакта 





[ато] т а ы М [ЕТ 











к К в и 
(12) 
1 
к - =" Ч&ат=8- А(х- с)? Ву?, (хе. 
После замен: 
Хх, =Х-С, Е, =&—с, 4. (х,,У) = 9(%,), 0, 0, 5, <> 5_ (13) 
ИУ (12) можно переписать в форме (звездочки далее опускаем): 
[|4 (Ельхь у) =в(т,у), (х,у)еО, (14) 
[®) 
1 1 х-& л+6+2с 
К я [2 — += = + 
($1,2) ов В? 
1 1 2 и 
-1т и у -Я + 1т[ х+ё+ р. 
А А А ^ й А 
&(х,у)=5- Ах’ Ву’, К =[(х+ё+2с)' +(у-п)"]. (16) 


При численном решении уравнения (14) применим метод нелинейных граничных ИУ типа 
Гаммерштейна, предложенный Галановым [7, 8]. Уравнение (14) дополним условиями неотрица- 
тельности контактного давления в области контакта, отсутствия контакта и обращения в нуль 
давления в дополнительной области 5\О, записав их все в виде системы 
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[ком м)а(мам = (М), 4(М)>0, М ео, 


[ку м)а(мам > (М), 4(М)=0, Ме$\О, —- 
5 
где введены обозначения М = (х, у), М = (&,"). 
Идея метода состоит в представлении искомого давления в форме 
4=9(М)=4'(М)+4 (М), (18) 
где введены нелинейные операторы 
4 (М) =зр{4"(М),0}, 4 (М)=шНа (М), 0}. (19) 


При учете (17) интегральное неравенство (16) будет удовлетворено в результате решения 
нелинейного операторного уравнения типа Гаммерштейна: 


Ор=0 (МЕС), Юр=р +К,р' -©, (20) 
где р. =р.(М), р; = р.(М), в=8(М), 
К.р’ = | К(М,М)р' (Мам. (21) 
5 


Можно доказать, что система (17) эквивалентна уравнению (20) [7]. Исследованы вопросы 
существования и единственности решения уравнения типа (20) [7]. Для численного решения 
уравнения (20) применим метод М. А. Красносельского, основанный на последовательных при- 
ближениях [7]. 

Прямоугольник $ покроем сеткой из л7 узлов с учетом отсутствия симметрии по координа- 
те х. Ясно, что с > е0, 50 < 1. Для проверки точности расчета ядра (15) можно использовать инте- 
грал [9] 

о р’ 
[ехр- Си), (Ри)аи = 


‚ П=О;1. 22 
о НЕЮ" (22) 


Исследуя поведение выражений Г, (и) —1 (п= 0;1) при и-> + для функций (9) - (11), по- 





ложим в (22) С= 2. При выполнении условия 


^> 41+ с) +1/(4 5) (23) 


для расчета ядра (15) можно использовать квадратурную формулу Гаусса по 32 узлам. 
Численные эксперименты. В табл. 1, 2 для задач А и Б соответственно даны значения кон- 
тактного давления 4(х,0) и вдавливающей штамп силы Р. Расчеты, результаты которых приведе- 
ны в таблицах, сделаны при у = 0,3, 6 = 0,004, А = 0,1, В = 0,005, = = с= 0,15 и разных значени- 
ях лин. 



































Таблица 1 
Значения давления и вдавливающей силы Р в задаче А 

9х0) х10? 3 
ь й х 0,075 —0,0375 0 0,0375 0,075 и 
1 0 1,77 3,00 3,39 3,32 2,81 0,640 
1 0,2 1,92 3,06 3,38 3,25 2,66 0,628 
1 0,2 1,64 2,94 3,39 3,39 2,95 0,651 
0,5 0 3,12 3,96 4,32 4,25 3,77 0,954 
0,5 0,2 3,23 4,00 4,31 4,18 3,64 0,941 
0,5 0,2 3,00 3,91 4,32 4,31 3,88 0,966 
0,3 0 4,63 5,42 5,70 5,55 4,99 1,42 
0,3 0,2 4,74 5,48 5,71 5,52 4,91 1,41 
0,3 0,2 4,51 5,35 5,67 5,58 5,07 1,42 
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Для задачи Б значения давления и силы меньше, чем для задачи А, как и должно быть. В 
обеих задачах контактное давление меньше на той стороне области контакта, которая ближе к 
участку, расположенному между штампами. 



































Таблица 2 
Значения давления и вдавливающей силы Р в задаче Б 

9(%0) х103 3 
ы н 0,075 —0,0375 0 0,0375 0,075 рн 
1 0 1,60 2,89 3,28 3,20 2,68 0,595 
1 0,2 1,74 2,95 3,27 3,13 2,53 0,585 
1 —0,2 1,46 2,82 3,27 3,27 2,82 0,605 
0,5 0 2,79 3,65 4,02 3,95 3,46 0,851 
0,5 0,2 2,92 3,70 4,01 3,89 3,34 0,842 
0,5 —0,2 2,66 3,59 4,01 4,00 3,57 0,858 
0,3 0 4,00 4,81 5,13 5,02 4,50 1,22 
0,3 0,2 4,13 4,89 5,17 5,02 4,46 1,23 
0,3 _0,2 3,88 4,73 5,07 5,01 4,53 1,21 
































Заключение. Решены новые пространственные контактные задачи с неизвестной областью кон- 
такта для упругого слоя о взаимодействии двух одинаковых эллиптических в плане штампов с 
учетом трения при различных граничных условиях на другой грани. При использовании метода 
нелинейных граничных интегральных уравнений с учетом симметрии задач определены области 
контакта, давления в этих областях, связи между силами и осадками штампов. Сделаны расчеты 
при разных значениях относительной ТОЛЩИНЫ СЛОЯ И коэффициента трения. 
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